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SUMMARY 

The Navy/NASA Engine Program (NNEP) Is a general purpose computer program 
currently employed by government. Industry and university personnel to simu- 
late the thermodynamic cycles of turbine engines. NNEP Is a modular program 
which has the ability to evaluate the performance of an arbitrary engine con- 
figuration defined by the user. 

In 1979, a program to calculate engine weight (WATE-2) was developed by 
Boeing's Military Division under NASA contract. This program uses a pre- 
liminary design approach to determine engine weights and dimensions. Because 
the thermodynamic and configuration information required by the weight code 
was available in NNEP, the weight code was appended to NNEP. 

Due to increased emphasis on fuel economy, a renewed interest has devel- 
oped in propellers. This report describes the modifications developed by NASA 
to both NNEP and WATE-2 to determine the performance, weight and dimensions of 
propellers and the corresponding gearbox. The propeller performance model has 
three options, two of which are based on propeller map interpolation. Propel- 
ler and gearbox weights are obtained from empirical equations which may easily 
be modified by the user. 


INTRODUCTION 

Due to the tremendous cost of turbine engine hardware, the ability to ana- 
lytically study the feasioiiity of advanced engine concepts is essential. One 
tool used to model the thermodynamic cycles of turbine engines is the 
Navy/NASA Engine Program (NNEP). NNEP has the ability to model arbitrary tur- 
bine cycles defined by the user. For design point calculations, the user 
inputs the engine configuration, component performance and optimization infor- 
mation. For jff-design calculations, component performance is determined 
using component maps. User-defined free variables are iterated to match com- 
ponent gas flows, rotational speeds, and power. For a further discussion of 
NNEP the user is referred to reference 1. 

In 1979, under NASA sponsorship, Boeing d(!ve1oped a computer program 
called WATE-2 that predicts weights of turbin^j engines (ref. 2). This program 
uses a preliminary design approach in which components are sized using geomet- 
ric and physical constraints along with thermodynamic cycle data. A logical 
source of this thermodynamic information is NNEP. Therefore, WATE-2 was writ- 
ten to be an integral part of NNEP but can be adapted to any thermodynamic 
cycle code. 



Renewed interest in propellers has grown rapidly due to high fuel prices. 
While NNEP is a very versatile program, its ability to model propellers is 
quite limited. The purpose of this paper is to describe modifications made to 
both NNEP and WATE-2 to determine propeller performance and weight. Three 
options of modeling propeller performance with varying degrees of sophistica- 
tion are available. The first option is simply a constant efficiency input. 
The second option interpolates from a generalized propeller map which is 
scaled for various integrated design lift coefficients, number of blades, 
activity factors and power coefficients. The third option interpolates from a 
user-supplied propeller performance map. 

Propeller weight and gearbox weight and dimensions are obtained from 
empirical relationships. These relationships are defaulted into WATE-2. The 
user may input his own empirical relations. The method is described later in 
the section on weight and dimension analysis. In the following discussions it 
is assumed that the reader has a basic understanding of NNEP and WATE-2. 
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SYMBOLS 


Blade activity factor 
Number of blades 

Power coefficient = 

p n D 

Thrust coefficient = 1 — j 

p D 
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Weight of propeller counterweights, lb 

Propeller diameter, ft 

Gear ratio = propeller speed/engina speed 

Propeller advance ratio = V/nD 

Airplane Mach number 

Propeller speed, rpm 

Propeller speed, rps 

Power loading, SHP/D^ 

Shaft horsepower 

Propeller thrust, lb 

Propeller tipspeed, ft /sec 

Flight velocity - true airspeed, ft/sec 

Component weight, lb 

Propeller efficiency, J(Ct^Cp) 

Mass density of air, slug/ft^ 
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The first option of propeller performance modeling uses a fixed propeller 
efficiency. The efficiency may be input as a different value from case to 
case but it is not dependent on propeller operating conditions. Propeller 
thrust is then simply calculated from 

T - 550 X SHP X n/V 

where n ii the input propeller efficiency and V is the forward velocity 
in tt/sec. 


For the static case this equation is indeterminate. Therefore, the thrust 
is calculated from the input static thrust-to-horsepower rat^o. 

T = (thrust/SHP) 5 ta^^Q x SHP 

Thus, only the shaft power, the efficiency, and static thrust-to-horsepower 
ratio are required fur this option. However, if the component weight is 
desired the propeller tip speed, power loading, number of blades, and activity 
factor are also required. 


The second option utilizes a generalized propeller map. This procedure 
was developed by Hamilton Standard under NASA contract (ref. 3). The advance 
ratio J and power coefficient Cp are first determined from the following 
equations: 


J = 


T.V 

T? 


where 


Cp 


550 SHP 

n3n5 
pn D 


D = /shpTdT 


TS 


Tnese values of Cp and J are used in a generalized map defined for the 
specified number of blades to obtain a thrust coefficient Ct. The general- 
ized maps have a base activity factor of 150 and a base integrated design lift 
coefficient of 0.5. Thus, the value of Ct is then modified through the use 
of adjustment factors to account for the effects of integrated design lift 
coefficient and activity factor. Finally, the propeller efficiency and thrust 
are determined by 



T = Ctpn^D^ 


For this option it is necessary for the user to know shaft horsepower, 
propeller power loading, number of blades. Integrated design lift coefficient, 
blade activity factor, and design tip speed. 
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The 9enera1ized maps were developed to be capable of predicting perform- 
ance for propellers with 2 to 8 blades, 80 to 20C blade activity factors, and 
0.3 to 0.8 integrated design lift coefficients. It should be noted that the 
original study only covered general aviation powerplants in the size range of 
100 to 1500 shaft horsepower. The ''alidity of performance predictions of pro- 
pellers used for powerplants outside this size range is unknown. 

The third option allows the user to input his own propeller map. Propel- 
ler maps are input with Ct as a function of Cp and J. The format of the 
map and a typical map are listed in appendix A. The format is the standard 
format required for all NNEP maps. The method of analysis is essentially the 
same for option 3 as it was for option 2. Advance ratio and power coefficient 
are first determined. These values are then used to obtain Ct from the 
input map. However, the value of Ct is not corrected for activity factor, 
number of blades, and integrated lift coefficient since these values are 
assumed implicitly fixed for a given map. Therefore, inputting these values 
has no effect on the performance calculated using option 3. The propeller 
effici( icy and thrust are then obtained the same as in option 2. 

In both option 2 and option 3 the user has two inputs available which 
allow the results predicted by the maps to be scaieu. One input allows the 
user to apply a constant scale factor to all thrusts and efficiencies pre- 
dicted by the maps. The other input allows the user to input a desired design 
point propeller efficiency. The program then sets up a scale factor which is 
the ratio of the desired design point efficiency to the map calculated effi- 
ciency. This scale factor is then applied to the map predicted performance in 
subsequent off-design cases. If both of these inputs are specified the con- 
stant scale factor overrides the desired design point efficiency. 

A list of the user-required inputs for each of the three options is given 
in table I. A more complete description of all the inputs is given in the 
User Manual (appendix B). 

While using any of the above options the user may desire a specified 
thrust from the engine. This may be obtained by the use of NNEP controls. To 
obtain the desired thrust for the design case the user may vary the inlet air- 
flow or the design turbine inlet temperature. For the off-design case the 
user may vary the turbine inlet temperature to obtain the desired thrust. 


PERFORMANCE PROGRAM STRUCTURE 

Propeller performance calculations are controlled by a subroutine called 
NNPROP. NNPROP is called from subroutine FLOCAL after the engine matching is 
complete. A copy of subroutine NNPROP is given in appendix C. 

When a constant efficiency is input (option 1) all propeller performance 
calculations are performed within NNPROP. Note that for this option NNPROP 
uses the same three inputs that were originally required by NNEP to run a pro- 
peVif.r. Therefore, any data sets which were previously used on NNEP should 
run with this new version without change. 

When the generalized maps are used (option 2) NNPROP makes a call to sub- 
routine PERFM. This subroutine contains the generalized maps and the correc- 
tion for integrated lift coefficient, number of blades, and activity factor. 
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The interpolation of the maps is done by two utility subroutines, BIQUAD and 
UNINT. A copy and a detailed description of these three subroutines are given 
in reference 4. 

For user-supplied maps (option 3), NNPROP calls the NNEP subroutine TLOOK 
to do the map interpolations. 


WEI3HT AND DIMENSION ANALYSIS 

WATE-2 was originally constructed to give weights and dimensions of compo- 
nents based on a preliminary design analysis. The routine to determine pro- 
peller and gearbox weight and dimensions, however, is based entirely on 
empirical data. 


Generalized propeller weight equations were obtained from Hamilton 
Standard, reference 5. One equation defines the weight of double-acting pro- 
pellers as follows: 



where Kw = 355. 


A second set of equations predicts the weight for single acting propellers 
with counterweights: 



To p*'ovide full flexibility, the propeller weight equation may be 
defaulted at the above values or the constant Kw and the parameter exponents 
may be input by the user. However, if the user desires to change any part of 
the equation, the constant and all the exponents must be input. It is also 
possible to. -nput a constant value for the counterweights C^, if desired. 


The value of Kw is defaulted at 355 for double-acting propellers and at 
220 for single-acting counterweighted propellers. These values imply current 
technology propellers using solid aluminum blades. The value of 355 repre- 
sents a higher speed application normally associated with double-acting pro- 
pellers, than the 220 value. Values of Kw ranging between 160 to 180 may be 
assumed for advanwed technology fiberglass or composite propellers. Other 
parameters required for the propeller weight equation, such as shaft hor:?- 
power, propeller rpm, etc., are obtained from maximum power conditions encoun- 
tered in NNEP. 
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There is also a provision for calculating gearbox weight. Figure 1 repre 
sents a defaulted curve used to estimate the gearbox weight as a function of 
output torque for a gear ratio of 0.118. For other gear ratios 

Wtgg = (Wtg0)unscaled 

The user may specify a different equation for weight as a function of 
torque by inputting the slope and y intercept. This curve must be input at 
a gear ratio of 0.118. 

Dimensions of the gearbox can also be calculated if desired. Defaulted 
values of these equations are 

Length = 0.125 (SHP)0-5 (1/GR)0.33^ in. 

Diameter = 0.230 (SHP)0*5 (1/GR)0.33, in. 

These equations are based on empirical data for in-line gearboxes in the 
1000 to 2500 horsepower range. All of the above equations may be modified by 
the user (see appendix B, User Manual). No provision is currently available 
for offset gearbox dimensions. 


WEIGHT PROGRAM STRUCTURE 

The control routine for WATE-2 is called WTEST. WTEST calls a subroutine 
(PROPWT) which performs the weight and dimension calculations. Copies of both 
routines are given in appendix C. Changes to the original WTEST are denoted 
by comment cards. 


CONCLUSION 

As interest increases in propeller research, analytical tools which accu- 
rately model both engine and propeller performance will be required. The 
modifications described herein to NNEP and WATE-2 provide such models. The 
program has been kept general enough to permit the user to update models as 
more detailed information becomes available. The weight and dimension calcu- 
lations, however, should be used with the understanding that the results are 
only rough estimates. It is suggested that the user examine these equations 
and make changes as required. 
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APPENDIX A 


PROPELLER MAP INPUTS 


Each propeller map has the following input data card setup. 

Card 1 Table reference number (Integer, col 2-5) Table identification label. 
Col 6-75 

Card 2 Z-Identifier (4 Character Symbol, col 1-4). Enter any dummy symbol 
because Z values are not used. 

NZ-Number of Z values (Integer, col 6 & 7). Must enter a value of I. 
Z-Variable values, 7F10., Beginning in Column 11. Must enter a 
value of 0. 

Card 3 Advance ratio identifier (4 Character Symbol, col. 1-4) 

NJ-Number of advance ratio values (Integer, col 6 & 7) 

XJ-Advance ratio values, 7F10., Beginning in Column 11. If needed, 
extra cards follow 10X,7F10. Format. 

Advance ratio values MUST be in ascending order. 

Card 4 Power coefficient identifier (4 Character Symbol, col 1-4) 

NCP-Number of Cp values (Integer, col 6 & 7) 

Cp values, 7F10., Beginning in Column 11. If needed, extra cards 
follow 10X,7F10. Format. 

Cp values MUST be in ascending order. 

Card 5 Thrust coefficient identifier (4 Character Symbol, col 1-4) 
NCP-Number of Cp values (Integer, col 6 & /) 

Ct values, 7F10., Beginning in Column 11. If needed, extra cards 
follow 10X,7F10. Format. These Values correspond to the values 
on the Cp Identifier Card. 

Last Card 3 Character Symbol EOT in col 1-3 

The propeller maps must be in the same dataset containing the other maps 
to be used by NNEP (i.e. turbine maps, compressor maps, etc.) 
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SAMPLE PROPELLER MAP 
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4100 H.S. ADV PROP 4 BLAOE 100 AF .55 CLI 


z 

1 

0.0 



J 

16 

0.0 

0.2 

0.4 



1.4 

1.6 

1.8 



2.8 

3.0 


CP 

10 

0.0568 

0.0737 

0.0926 



0.2747 

0.3192 

0.3558 

CT 

10 

0.1456 

0.1732 

0.1965 



0.2559 

0.2565 

0.2541 

CP 

10 

0.0499 

0.0655 

0.0841 



0.2406 

0.2818 

0.3200 

CT 

10 

0.1084 

0.1354 

0.1626 



0.2580 

0.2605 

0.2600 

CP 

11 

0.0406 

0.0564 

0.0769 



0.2234 

0.2673 

0.3840 

CT 

11 

0.0622 

0.0931 

0.1237 



0.2371 

0.2470 

0.2620 

CP 

10 

0.0369 

0.0564 

0.0819 



0.2561 

0.3007 

0.3200 

CT 

10 

0.0362 

0.0712 

0.1051 



0.2439 

0.2551 

0.2600 

CP 

9 

0.0482 

0.0797 

0.1156 



0.3347 

0.3857 


CT 

9 

0.0447 

0.0828 

0.1186 



0.2572 

0.2649 


CP 

10 

0.0323 

0.0697 

0.1141 



0.3722 

0.4264 

0.4821 

CT 

10 

0.0158 

0.0585 

0.0985 



0.2561 

0.2694 

0.2696 

CP 

8 

0.0536 

0.5345 

0.1070 

0.1639 

CT 

8 

0.0341 

0.2782 

0.0785 

0.1194 

CP 

9 

0.0337 

0.0965 

0.1654 



0.5307 

0.5997 


CT 

9 

0.0108 

0.0600 

0.1052 



0.2754 

0.2883 


CP 

9 

0.0695 

0.1337 

0.2011 



0.5395 

0.6084 


CT 

9 

0.0339 

0.0744 

0.1121 



0.2624 

0.2823 


CP 

8 

0.0407 

0.5815 

0.1134 

0.1913 

CT 

8 

0.0094 

0.2588 

0.0543 

0.0957 

CP 

6 

0.0937 

0.1836 

0.2749 

CT 

6 

0.0368 

0.0821 

0.1238 

CP 

6 

0.0773 

0.1797 

0.2853 


0.6 

0.8 

1.0 

1.2 

2.0 

2.2 

2.4 

2.6 

0.1170 

0.1484 

0.1866 

0.2287 

0.2179 

0.2327 

0.2456 

0.2531 

0.1058 

0.1299 

0.1614 

0.1999 

0.2890 

0.2083 

0.2300 

0.2480 

0.1001 

0.1253 

0.1547 

0.1885 

0.3200 

0.1532 

0.1812 

0.2075 

0.2290 

0.2640 

0.1110 

0.1421 

0.1758 

0.2133 

0.1376 

0.1685 

0.1972 

0.2229 

0.1541 

0.1940 

0.2378 

0.2850 

0.1525 

0.1844 

0.2136 

0.2389 

0.1607 

0.2093 

0.2597 

0.3149 

0.1355 

0.1706 

0.2025 

0.2321 

0.2222 

0.2821 

0.4117 

0.4777 

0.1527 

0.1930 

0.2533 

0.2740 

0.2346 

0.3055 

0.3768 

0.4542 

0.1462 

0.1850 

0.2188 

0.2505 

0.2657 

0.3323 

0.4006 

0.4667 

0.1466 

0.1796 

0.2105 

0.2364 

0.2692 

0.3464 

0.4246 

0.5015 

0.1338 

0.1691 

0.2026 

0.2318 

0.3643 

0.4542 

0.5855 


0.1617 

0.19/6 

0.2308 


0.3891 

0.4916 

0.5951 
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CT 

6 

0.0227 

CP 

7 

0.0675 

CT 

7 

0.0130 

CP 

7 

0.0968 

CT 

7 

0.0223 

CP 

8 

0.0759 

0.6046 

CT 

8 

0.0080 

0.1909 

CP 

7 

0.0996 

CT 

EOT 

7 

0.0132 


0.0720 

0.1170 

0.1578 

0.1531 

0.2433 

0.3341 

0.0532 

0.0905 

0.1254 

0.1794 

0.2649 

0.3510 

0.0576 

0.0904 

0.1213 

0.1473 

0.2230 

0.3005 

0.0387 

0.0674 

0.0947 

0.1810 

0.2666 

0.3537 

0.0451 

0.0754 

0.1039 


0.1958 

0.2312 


0.4231 

0.5103 

0.5980 

0.1576 

0.1878 

0.2167 

0.4361 

0.5188 

0.6010 

0.1503 

0.1773 

U.2031 

0.3784 

0.4556 

0.5312 

0.1209 

0.1457 

0.1689 

0.4410 

0.5272 

0.6112 

0.1309 

0.1565 

0.1804 
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APPENDIX B 


USER MANUAL 

The first part of this manual describes the inputs required to run propeller 
performance calculations, while the later part presents the weight code 
inputs. It is assumed that the user is familiar with NNEP and its method of 
input. For instructions the reader is directed to "K0NFIG/REK0NFI6" (ref. A) 
which is an interactive input preprocessing program. A sample input d..t 
is given at the end of this Appendix. 

The NNEP input string used to read in propeller performance data is the 
KONFIG statment. All propellers will be read in the following form: 

K0NFIG(1,N)='L0A0',SPEC(1,N)= (1),(2) (10) 

where N is the component number. 
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INPUT PROPELLER ARRAY 


Spec(l ,N) 
Spec(2,N) 


Spec(3,N) 

Spec(4,N) 

Spec(5,N) 

Spec(6,N) 

Spec(7,N) 

Spec(8,N) 


Spec{9,N) 


Shaft horsepower (negative) delivered to the propeller 
One of three options: 

1) constant efficiency - if a constant efficiency is specified, 
the value must be 0 < n < 1 

2) negative value - any input, negative in sign, will initiate 

the generalized propeller table interpolation (PERFM). 
Example: -.1, -100, etc. Input magnitude has no 

significance. 

3) Table number - table number of user input propeller 
performance table. This value must be <1. 

Note: loads other than propellers may be connected to shafts 

in NNEP engine configurations; for such loads SPEC(2) must be 
set to 0.0. 

Propeller thrust divided by shaft horsepower (T/SHP) for stacic 
thrust calculations.^ 

Design point power loading, PL = SHP/D^.b 
Number of blades of the propel ler.c 
Integrated design lift coefficient.*^ 

Blade activity f actor. c 

A constant scaling factor which is multiplied by the results 
taken from the propeller map. This input scales the efficiency, 
thrust, and coefficient of thrust when using option (2) or 
option (3) of Spec(2,N). Spec(8,N) will override Spec(9,N) 
if both are input non-zero. 

Desired propeller design point efficiency. This input also 
determines a scale factor which will be applied to all 
subsequent off -design cases. 

Spec(8,N) will override Spec(9,N) if both are input non-zero. 


Spec(10,N) =Design tip speed of propeller (ft/sec). ^ 


^This input is only used for option(l) of Spec(2,N) 

t>This input used to get prop diameter and rpm which are required for options 
(2) and (3) of Spec(2,N) and for propeller weight calculations. 

^This input only used for option (2) of spec(2,N) and for propeller weight 
calculations. 

<^This input only used for option (2) of spec(2,N). 
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OUTPUT PROPELLER ARRAY 


Datout(l,N) 

Datout(2,N) 

Datout(3,(l) 

DatouH4,N) 

0atout{5,N) 

Datout(6,N) 

Datout(7,N) 

Datout(8,N) 

0atout(9,N) 


Shaft horsepower delivered to the propeller. 

■’ropeller rotational speed (rpm). 

Propeller thrust, (Ibf). 

Advance ratio, J. 

Coefficient of power, Cp. 

Coefficient of thrust, Ct. 

Propeller tip speeo, (ft/sec). 

Unadjustt^ p»'opfc'iler efficiency. 

Adjusted propeller efficiency, 

(using SPEC(8,N) or SPEC(9,N)). 

MECHANICAL DESIGN INDICATORS 


The mechanical design indicators (IWMEC) are required to determine whether the 
propeller and gearbox are to be weighed. IWMEC is a two dimensional integer 
array, and is of the form IWMEC(M,N) where N is the component number of the 
prop as used in NNEP, and M is the variable identifier as defined below. 

For a propeller: 


IWMEC ARRAY 

location M DESCRIPTION 


1 Type of component being weighted. 

'PROP' -- indicates propeller. 

2 Indicates whether propeller weight is to be calculated. 

0 -- Do not calculate propeller weight. 

1 -- Calculate propeller weight. 

1 Indicates presence of counterweights on propeller. 

0 -- No counterweights (double acting prop). 

1 -- ProD contains counterweights. 

4 Indicates whether gearbox weight is to be calculated. 

0 -- Do not calculate gearbox weight. 

1 -- Calculate gearbox weight. 

5 Indicates whether gearbox dimensions are to be 
calculated. 

0 -- Do not calculate gearbox dimerisions. 

1 — Calculate gearbox dimensions. 

6 Blank 


7 Blank 
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Notes on IWMEC Inputs 


If IWMEC(2,N) is input as 0, a message is printed informing the user that 
propeller weight will not be calculated. If input is 1, the exponents of the 
propeller weight equation are alterable. See OESVAL array discussion. 

If IWMEC(4,N) is input as 0, a message is printed informing the user that 
gearbox weight will not be calculated. 

If IWMEC(5,N) is input as 0, a message is printed informing the user that 
gearbox dimensions will not be calculated. It is assumed that the propeller 
adds no length to the engine beyond that of the gearbox. 



ILENG Array 


The ILENG array indicates which components will be included in the length 
summation. The propeller component number should be included to add the 
length of the gearbox to the engine. 

DESIGN VALUES 


This section describes mechanical and aerodynamic design data necessary to 
determine the weight and dimensions of the propeller. The data required is in 
the floating-point two-dimensional array DESVAL(M,N), where M is the variable 
array location identifier. N is the component number of the prop as used in 
NNEP. 


The following equations are used to calculate propeller weight and gearbox 
weight, length, and diameter: 

fn OESVAL(IO) R DESVAL(ll) .c DESVAL(12) 

'^^PROP DESVAL(2) ^(^0) (q-) (’y^) 


.. DESVAL(14) 5.p DESVALdS)] 

ND ^ /nccwAi ni4.11 J + 


DESVAL(13) 

W (DESVAL(l).l) 


10D‘ 


Cw 


Cw = 5 




^ . 20000 .^*^ 

^ ND ’ 


1 


or if DESVAL(7) > 0.0 
Cw = DESVAL(7) 


Wtgg = DESVAL{4) * [dESVAL( 5) * TORQUE + DESVAL(6)j 


18 

W 


DESVAL(17) 


) 


(IB) 

(2B) 


(3B) 


LENGTHgg = DESVAL(8) (SHP)°*^ {1/GR)°‘^^ (^B) 

DIAMETERgg = DESVAL(9) (SHP)^‘^( 1/GR)°‘^^ (5B) 
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The DESVAL variables used in the above equations are described as follows: 
DESVAL ARRAY 

LOCATION DESCRIPTION 

(1) Design Mach number of aircraft (see equation IB). Default input of 
zero yeilds zero Mach number input. 

(2) Kw - Coefficient of propeller weight equation (see equation (IB)). 

A multiplier representing the level of propeller technology, 
described in the body of this report. If Kw is input as non-zero, 
the exponents of the weight equation (DESVAL(10,N)-(15,N)) must be 
input. If Kw is defaulted (zero), default values are used. 

(3) Gear ratio (prop rpm/shaft rpm); see equations 3B,4B, and 5B. 

OR Shaft RPM. The user can choose whether to specify a gear ratio 
(usually < 1) or an RPM of the shaft delivering power to the prop. 

This RPM is independent of the RPM calculated by NNEP. If 
0ESVAL(3,N) > 200. The input is assumed to be RPM and the 
corresponding gear ratio is calculated and stored in DESVAL(3,N). 

If DESVAL(3,N) < 200. The input is assumed to be the gear ratio. 

(4) Coefficient of gearbox weight equation (see equation 3B). No scaling 
will occur if 0.0 or 1.0 is input. 

(5-6) Slope and Y-Intercept, respectively, of optional user input 
gearbox weight curve (see equation 3B). If slope is input a 
Y-Intercept must also be input. Any input curve must be scaled at 
0.118 gear ratio. To use default curve described in the body of 
the report input 0.0 for both values. 

(7) Input value of counterweight weight (see equation 2B). The default 
calculation value may be used by inputting DESVAL(7,N)=0.0. 

(8) Scale factor for gearbox length (see equation 4B). No scaling will 
occur if 0.0 or 1.0 is input. 

(9) Scale factor for gearbox diameter (see equation 5B). No scaling will 
occur if 0.0 or 1.0 is input. 

(10-15) Exponents of propeller weight equation (see equation IB). The 

default values are used if 0ESVAL(2,N) = 0.0, If DESVAL(2,N) =0, all 
of the exponent values must be input. 

(16) Not used. 

(17) Exponent on gear ratio term of gearbox weight equation (see 
equation 38). 
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DESVAL DEFAULT VALUES 


DESVAL VALUE ABBREVIATED 
DEFAULT ARRAY DESCRIPTION 
LOCATION 


1 

0.0 

Design Mach numUer 

2 

355.0 

(220.0)3 

Kw— prop weight scalar 

3 

1.0 

Gear ratio 

4 

1 .0 

Gearbox weight scalar 

5 

0.0174 

Gearbox weight curve slope 

6 

45.0 

Gearbox weight curve 
Y-Intercept 

7 

0.0 

(equ. 2B)3 

Weight of counterweights 

8 

1.0 

Gearbox length scalar 

9 

1 .0 

Gearbox diameter scalar 

10 

2.0 » 


11 

0.7 


12 

0.75 j 



(0.70)3 

Prop weight equation 
> exponents 

13 

0.5 ( 

(0.4)3 1 


14 

0.5 


15 

0.12 ; 


16 

— 

Not used 

17 

0.5 

Exponent on gear ratio 
term in equation 3B. 


^Defaults for propellers with counterweights. 
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SAMPLE INPUT CASE 


ORIGINAL PAGE IS 
OF POOR QUALITY 


THIS IS A TITLE CARD 
&0 CALGLO«T,AMAC=r &cND 
iO M00E=1,JWT=1 

K0NFIG(1.1)=4HINLT,1,0,?.0,SPEC(1,1)=14. 3.0, 0,0,. 23, 0,0, 0,5000, 
K0NFIG(1,2)=4HC0MP, 2,0, 5.0, SPEC(1,2)=1. 278, 0,1, 3707, 1,3708,. 982, 3709, 1,0,0,. 87, 

5.896. . 919, 

KONFIG( 1,3)=4HC0MP, 3, 0,4, 20, SPEC(1,3)«1. 803,. 05, 1,3404, 1,3405,. 78. 3406, 1,0,0, 

.864,3.887,0.988, 

KONFIGI 1,4)=4HDUCT, 4, 0,5,0, SPEC! 1.4)«. 05,0, 0,2710,. 985, 18300, 

KONFIG( 1,5)=4HTURB, 5, 20, 6,0, SPEC(1,5)=3. 5,. 95, 1,3801, 1,3802,. 8881,1,. 85,1,. 91, 

5665,1,08,1, 

KONFIG( 1,6)=4HTURB, 6, 20, 7.0, SPEC(1,6)=2. 49,. 05, 1,3803, 1,3804,. 9154, 1,1,1,. 9156, 

5124.1.000000. 3, 

KONFIG( 1, 7)=4HN0ZZ, 7, 0,8, 0,SPEC(1.7)=30, 1,0,0,. 95, 3*0,1, 

K0NFIG(1,8)»4HL0AD,SPEC(1,8)=-2686,.985.3.5,9.6,4,.5,180,1.0115,0.00,900., 

K0NFIG(1,9)=4HSHFT,2,3,5,0,SPEC(1.9)=38834,4*1,2*1,.99, 

K0NFIG(1,10)=4HSHFT,8,11,6,0,SPEC(1,10)»18000,.0611,1,1,1..97,1,1 

K0NFIG(1,15 =4HCNTL,SPCNTL(1.15)=1,1,4HSTAP,8,2,0 

KONFIG( 1,11 )=4HL0AD,SPEC( 1 ,1 1 )=-68, 

KONF IG ( 1,16 =4HCNTL , SPCNTL ( 1 . 1 6 ) =1 , 2 .4HSTAP ,8 , 3 , 0 
K0NFIG(1,17)=4HCNTL,SPCNTL(1,17)=1,3,4HSTAP,8,5,0 
K0NFIG{1,18)=4HCNTL,SPCNTL{1.18)*1,5,4HSTAP,8,6,0 
KONF I G ( 1 , 1 9 ) =4HCNTL , SPCNTL ( 1 , 1 9 ) = 1 ,6 .4HSTAP ,8,7,0 
K0NFIG(1,20 =4HCNTL,SPCNTL(1,20 =1,9,4H00UT,8,9,0 
KONFIG( 1 ,21 )=4HCNTL,SPCNTL( 1 .21 )=1 ,8,4HD0UT,8, 10,0 
KONFIG(l,22 =4HCNTL, SPCNTL! 1,22 =2,3, 4HPERF, 15, 0,0,1 
K0NFIG(1,23)=4HCNTL.SPCNTL(1,23)=10,2,4HD0UT,5,2,15,0,-5.10, 

&ENO 

40 CALBLD=F,SPEC(9,22)=0,SPEC(9,15)»1,SPEC(9,16)=1,SPEC(9,17)=1,SPEC(9.18)=1, 

IWT=1,SPEC(5,1)=0,SPEC(9,1)=0,SPEC(4,4)=2760,SPEC(9,19)«1,SPEC(9,20)=1. 

SPEC(9,21)=1, &ENO 

40 SPEC(5,1)=0.0,IWT=1, 4EN0 

40 SPEC(5,1)=.1 .IWT=2, 4EN0 

4W ISII=F,ISI0=F,I0UTC0=2,0ISKWI=1 

1LENG=8,2,3,4,5.6,7, 

IWMEC(1,2 ='LPC',1,1,0,0,0,3 
IWMEC(1,3)='HPC',2,0,1,0,0,1 
IWMEC 1,4)='PBUR',1 
1WMEC(1,5)='HPT', 1,2, 0,1,0 
IWMEC(1,6)='LPT', 1,2, 0,3,0 
IWMEC(1,7)»'N0ZZ', 1,0,0 
IWMEC(1,8 *'PR0P',1,0,1,1 
1WMEC(1,9)='SHFT',2,5,0,0,2 
IWMEC(1,10)='SHFT',1,6,0,0,8 

DESVAL(1,2)«. 55, 1.9, 0.5, 1.5,4. 0,3. 0,0. 45,0.0,0.0,1.0,0.0.1.0,1.0,4*0.0, 

DESVAL! 1,3 =0.4, 1.4, 0.7, 1.5, 3. 0,1. 5,0. 3, 0.0, 0.0, 1.0, 0.0, 1.0, 1.0, 4*0.0, 

DESVAl (1.4) = 150., 0.015, 13*0.0 
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1)ESVAL{1 

nhSVAL(l 

UtSVAL(l 

UtSVAL(l 

L'ESVALfl 

OESVAL(i 

&END 


.U))=60(XW..0.^86,0.86,M;i^lo?6 
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APPENDIX C 


ORIGINAL PAGE (S 
OF POOR QUALITY 


PROGRAM LISTING 


c 

C SUBROUTINE NNPROP 

C 

C PURPOSE: 

C 
C 

C DATE OPERATIONAL: 

C 

C AUTHOR: 

C 

C DESCRIPTION: 

C 
C 
C 
C 
C 

(^**A****************************************1^********************************** 
^^****************************************************r************* ************* 

SUBROUTINE NNPROP (VKTS.ALT) 

IMPLICIT REALMS (A-H.O-Z) 

DIMENSION ZMS(2) 

COMMON /OBL / DAT I NP ( 1 5 , 60 ) , DATOUT ( 9 , 60 ) , WTF ( 40 ) , TOPRES ( 40 ) , TOTEMP ( - 
140),FAR(40),C0RFL0(40),VMACH(40).STATP(40),ERR0R(40),T0L,T0LT,T0LT- 
2 T , DE PV ( 20 ) , DTOL ( 20 ) , PERPF ( 20 ) , RCH , STOC , TFUE L 

COMMON /SNGL/ JMl , JM2 , JP 1 , JP2 , JCX, LOCTBL ( 9 , 60 ) , JCOMP ( 70 ) , I WAY , N I T, - 
^TAB(70), JC0NF(60,4), jrYPE(60), JFL0W(70),IDEDAP( 15),KKINDS(14,25),- 
2NC0MP , NOSTAT, N I TER , NF I N I S. NPASS ,JCC,NTBL,NCTS,JCIND(20),JCDEP(20),- 
3 JC V I ND ( 20 ) , JC VOE P ( 20 ) , KDT YP ( 20 ) , I DONE ( 60 ) 

N11=KKINDS(11,1) 

DO 2 IK=1,NH 
JXSHFT=KKINDS(11,IK+1 ) 

1)0 2 IISHFT=1,4 

IF(JCX.EQ.JCONF(vlXSHFT,IISHFT))GO TO 5 


THIS ROUTINE CALCULATES PROPELLER PERFORMANCE 
USING ONE OF THREE OPTIONS 

08/13/82 

R. M. PLENCNER 


SUBPROGRAM CALCUALTES PROPELLER THRUST FROM 
ONE OF THREE OPTIONS: 

1 INPUT CONSTANT EFFICIENCY 

2) GENERALIZED HAMILTON STANDARD PROPELLER MAP 

3) USER SUPPLIED PROPELLER MAP 
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2 CONTINUE 

5 DATINP(15,JCX) = JXSHFT 
PTRPM*0AT0UT(2,JXSHFT) 
SHP = -DAT0UT(1,JCX) 

PL = DATINP(4,JCX) 

BLAOT = DATINP(5,JCX) 
CLI = DATINP(6,JCX) 

ATF = DATINP(7,JCX) 

XFT = 0ATINP(8,JCX) 
ETADES = DATINP(9,JCX) 
TSDES = DATINP(10,JCX) 
LTAB=0 


ORIGINAL PAGE IS 
OF POOR QUALITY 


C THIS IS TO IMPLEMENT THE PERFM SUBROUTINE W/ GENERAL MAPS— 
IF(DATINP(2,JCX).LT.O.O)LTAB*-777777. 

C — IF TABLE NO. IS SPECIFIED, SKIP NEXT 4 LINES — 

IF(0ATINP(2, OCX). LE. 0.0. OR. DATINP(2,JCX).GT. 1.0)60 TO 633 
ETA1*0ATINP(2,JCX) 

IF(VKTS.EQ.0.0)ETA1=0.0 

C — SET ETAl EQUAL TO SPECIFIED ETA.CALCULATE THR FOR STATIC CASE — 
IF(VKTS.EQ.0.0)THR=DATINP(3,JCX)*SHP 
633 IF(LOCTBL(2,JCX).GT.O) LTAB=L0CTBL(2, JCX) 

C — 60 TO 10 TO IMPLEMENT SUB. PERFM — 

IF(LTA8.NE.0.0)G0 TO 10 

C — IF DATINP(9) DOES NOT AGREE W/ ETAl, WRITE AN ERROR MESSAGE — 

IF(NIT. LE.l. AND. ETADES. GT. 0.0. AND. ETAl. NE.ETAUES)WRITE(10, 300) JCX 

C — CALCULATE THRUST FOR NON-STATIC CASE — 

IF(VKTS.GT.O.O) THR=550*SHP*ETAl/( VKTS*1 .6878) 

THR1=THR 


IF(VKTS.GT.0.0)ETA=ETA1 
IF(PL.EO. 0.0. OR. TSDES. EQ.0.0)G0 TO 101 
C — CLACULATIONS BELOW WILL NOT BE REQUIRED FOR SPECIFIED 

C ETA OR TABLE LOOK-UP CASES 

10 CALL ICAO (ALT, TP, PRESS, SSFPS.DNSTY) 

DRAT = .002377/DNSTY 
IF (IWAY.EQ.l) DIA = DSQRT(SHP/PL) 

IF (IWAY.EQ.l) TS=TSDES 

IF (IWAY.EQ.l) DESRPM=60.0*TS/(3.14159*DIA) 

IF (IWAY.EQ.l! RPMRAT=OESRPM/PTRPM 

RPM=RPMRAT*PTRPM 

TS = 3.14159*RPM*DIA/60. 

IW = 1 

CP = SHP*10.E10 * DRAT/(2. * TS**3 * DIA**2 * 6969.26) 

ZJI = 101.27 * VKTS/(RPM * DIA) 

ZMS(l) = VKTS* 1.6878/SSFPS 
20 ZMS(2) = TS/SSFPS 
ZZ=0.0 

IF(LT<\B.LT.0.0) CAL! PERFM (IW,CP,ZJI,AFT,BLADT,CLI,CT1 ,ZMS,LIMIT) 
IF(LTA'B.GT.0.0 CALL TLOOK (LTAB,CP,ZJI,ZZ,CT1 ) 

IF(LTAB.NE.0.0)THR1 « CTl* TS**2 * DIA**2 * 2.4084E-4/DRAT 
IF(LTAB.NE.0.0)ETA1 = ZJI*CT1/CP 

IF(NIT. LE.l. AND. ETADES. NE. 0.0. AND. ETAl. EQ. 0.0) WRITE(10,100) JCX 
IF(NIT. LE.l .AND. ETADES. NE.0.0. AND. XFT. NE.0.0. AND. XFT. NE.l .0)- 
1WRITE(10,200)JCX 
IF(IWAY.NE.1)60 TO 40 
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ORIGINAL PACE IS 

ETAFIX=ETADES OF POOR QUALITY 

I F ( ETADES . EQ . 0 . 0 ) ETAF I X=ETA1 
ETARAT=1 .0 

I F ( ETAl . NE . 0 . 0 ) ETARA7 =ETAF I X/ETAl 
IF(XFT.NE.1.0.AN0.XFT.NE.0.0)ETARAT=XFT 
40 ETA=ETARAT*ETA1 
THR = THR1*ETARAT 

CT = THR*DRAT/(TS**2 * DIA**2 * 2.4084E-4) 

101 DAT I NP ( I I SHFT+1 . JXSHFT ) =RPM/DATOUT ( 2 , JXSHFT ) 
DAT0UT{IISHFT+2,JXSHFT)=RPM 
0ATINP(3,JCX) = THR/SHP 
DAT0UT(2,JCX) = RPM 
DAT0UT(3,JCX) = THR 
DAT0UT(4.JCX) = ZJI 
DAT0UT(5,JCX) = CP 
DAT0UT(6,JCX) = CT 
DAT0UT(7,JCX) = TS 
DATCUT(8,JCX) = ETAl 
DAT0UT(9,JCX) = ETA 

lOOFORMATC ','*** WARNING *** FOR L0AD#M2,' 0ATINP(9) I6N0RED- 
1 BECAUSE VEL0CITY=0 FOR DESIGN CASE') 

200 FORMATC ',10X,'*** WARNING ***'/' BOTH DATINP(8) AND (9) HAVE- 

1 BEEN SPECIFIED FOR LOAD #',I2,/' DATINP(8) WILL OVERRIDE DATINP(9)') 
300 FORMATC ','*** WARNING *** FOR LOAD#', 12,- 

1' DATINP(9) NOT = DATINP(2) -| DATINP(2) OVERRIDES ') 

RETURN 

END 


★ ★★★ 

C 

C SUBROUTINE PROPWT 

C 

C PURPOSE: 

C 
C 

C DATE OPERATIONAL: 

C 

C AUTHOR: 

C 

C DESCRIPTION: 

C 
C 
C 
C 
C 


c********************************************************* 


THIS ROUTINE CALCULATES THE WEIGHT OF A PROPELLER 
AS WELL AS WEIGHT OF AN APPROPRIATE GEARBOX. 

08/13/82 

PETER SENTY.NASA SUMMER EMPLOYEE 1982 

SUBPROGRAM UTILIZES GEOMETRICAL DATA FROM 
BOEING "WATE-2" PROGRAM ASWELL AS THERMODYNAMIC 
DATA FROM THE KAVY/NASA ENGINE PROGRAM TO PREDICT 
THE WEIGHT OF A PROPELLER, WEIGHT OF A COMPATIBLE 
GEARBOX AND THE DIMENSIONS OF THAT GEARBOX. 


it'k'kifif'k'kif'kidf'kitifififif'k'kifif'kirkirk'kif'kif'kit 


SUBROUTINE PROPWT (NC) 

REALMS DAT I NP , DATOUT , WTF , TORRES , TOTEMP , FAR , CORFLO , VMACH , STATP , ERRO- 
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o o o <-> o cn C-> 


c 

c 


IP TOi Tni T TOI TT DEPV QT0L.PERPF,RPMNT,TMTEMP»TMPRES,DAT0UM,DATMAC“ 

^bKokLIMJNPRES.^^^ 

3 ! DEBUG , OEPQ , SELAST , DD , TOLOPT 
***************** 

* COMMON BLOCKS * 

***************** ,/vTf-un/ 

COMMON /DBL/ DATINP(15,60) ,DAT0UT(9,60) ,WTF(40) J0PRES(40) 

1 40 ) , F AR ( 40 ) , CORFLO ( 40 ) , VMACH ( 40 ) , STATF ( 40 ) , ERROR ( 40 ) , TOL , TOLT , TOLT 

2UTCD,NSTAG(60) 

SS /NEPOPT/ SKdEpIsELAST.DDJOLOPT.NDSET.NPAR^ 
^C0MM0N^/TERM0N/*TNPRES(40) JNTEMP(40) ,CNRFL0(40) ,WN^^ 

l),DATOUN(9,6O).FARN(4O),0ANINP(15,6O),TMTEMP(4O),WRES{4O),CORFLM- 

2(40),DAT0UM(9,.60),DATMAC(4,60),DATALT(4,60),DESLIM(15) 

COMMON /CENTER/ CGARM(60) 

COMMON /ARM/ ACCARM 
COMMON /DISKK/ DISKWI, ENGINE 
COMMON /SKIP/ SKIPIT 

*************************** 

* DATA STORAGE DEFINITION * 

*************************** 

LOGICAL PINP,IPLT,ISIO,ISI I, PLOT, SKIPIT 
INTEGER IDID(60),ILENG(40) rr.ar^riAn\ 

DIMENSION NUMNUM(17), IRNAME(17), C0RFLC(40) Tm-rm ti fnp nF«; 

NAMELIST / W/IWMEC,DESVAL,ACCS,IWT,IPLT,ISII,ISIO,IOOTCD,ILENG,DES- 

1 L IM , I SCALE , SCALE , ACCARM , 0 I SKW I , ENG I NE , PLOT 
DIMENSION X(17), DEFALT(17) 

★*★*★★*★★★*★★★★★*** 

* DATA STATEMENTS * 

******************* 

DATA I DUC , LSHAF , ENGU , S I U , I VAL V/4HDUCT .4HSHAF * ^^HF^SFy- 

DATA ILPC,IHPC,IFAN,IF0.IFI,IHPT,ILPT/3HLPC,3HHPC.3HFAN.2HF0.?HFI. 
^DATA^IPR0P/4HPR0P/,DEFALT/0.0,355.,1.0,1.0,0.0174,45.0,0.0,1.0,- 

11 .0,2. 0,0. 7,0. 75,0. 5, 0.5, 0.12, 220., 0.5/ 

C TEST CALL TO PROPWT — 

IF(IWMEC(1,NC).NE.IPR0P) RETURN 
C — INITIALIZE "X'' ARRAY- 
00 5 1=1,17 
X(I)’0.0 
5 continue 

C— FILL X ARRAY WITH INPUTS OR DEFAULT VALS.— 

IF{0ESVA^O,NC).NE.0.0) X( I )*DESVAL( I .NC) 

IF(DESVAL(I,NC).EQ.O.O) X( I )»DEFALT( I ) 

10 CONTINUE 


ORIGINAL PAGE IS 

OF POOR QUALITY 
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OR’WP..’. JS 
or POOR QO-vLITY 
C 

C ******* BEGIN CALC. OF GEAR BOX WT. & DIMENSIONS ***** 

C — attain max. value PARAMETERS — 

SHP=DATOUM(l,NC) 

RPMAX=DAT0UM(2,NC) 

C --- LOCATE SHAFT — 

JSHP=0ATINP(15,NC)+.01 
C — CHECK FOR GEAR RATIO OR RPM INPUT — 

IF(DESVAL(3,NC).NE.0.0.AN0.0ESVAL(3,NC).LT.200.) S6R=DESVAL(3,NC) 
IF(DESVAL(3.NC).NE.0.0.AND.DESVAL(3,NC).GT.200.)- 
1 SGR=RPMAX/DESVAL(3,NC) 

IF(DESVAL(3,NC).NE.0.0)G0 to 35 

RSMAX=RPMAX/SGR 

DO 500 IX=1,4 

I F { NC . EQ . JCONF ( JSHP , I X ) ) SGR=DATOUT ( 2 , JSHP ) 
WRITE(10,110)SGR,RPMAX,SHP 
500 CONTINUE 
35 CONTINUE 

C — SEND RPM BACK TO SUB. "TURBINE" FOR USE IN DIMENSION CALCS — 
JSHP=DAT0UT(15,NC) 

RPMT(NC)=RSMAX 
CALCULATE TORQUE — - 
T0RK=SHP*5252./RPMAX 
C — CHECK OPTION TO WEIGH G.B. - — 

IF(IWMEC(4,NC).EQ.O) WRITE(10,210) 

IF{IWMEC(4,NC).EQ.O) GO TO 100 
C — INITIALIZE VALUES TO BE PRINTED OUT — 

GBWT=0.0 

GBKW=DEFALT(4) 

I F { OES VAL ( 4 , NC ) . NE . 0 . 0 ) GBKW=DESVAL ( 4 , NC ) 

C 

C — CHECK OPTION FOR G.B. WT. CURVE — 

X(5)=0ESVAL(5,NC) 

X(6)=0ESVAL(6,NC) 

IF(DESVAL(5,NC).EQ.0.0.AND.DESVAL(6,NC).EQ.0.0) X(5)=0EFALT(5) 
IF(0ESVAL{5,NC).EQ.Q-C.AN0.DESVAL(6,NC).EQ.0.0) X(6)=0EFALT(6) 
GBM=X{5) 

GBB=X(6) 

C — CHECK FOR INPUT GEAR RATIO SCALAR EXPONENT — 
IF(0ESVAL(17,NC).NE.0.0) X(17)=DESVAL(17,NC) 
IF(OESVAL(17,NC).EQ.O.O) X(17)=0EFALT(17) 

C— INTERPOLATE CURVE FOR GEARBOX WEIGHT — - 

GBWT=G8KW* ( GBM*T0RK+C3B )*((!. /SGR ) /8 . 5 ) **X( 1 7 ) 

100 CONTINUE 

C — CHECK FOR INPUT SCALE FACTORS FOR G.B. DIMENSIONS — 

GBLN=0.0 

GB0I=0.0 

IF(IWMEC(5,NC).EQ.0) WRITE(10,220) 

IF(IWMEC(5,NC).EQ.O) GO TO 101 
X(8)=0EFALT(8) 

X(9)=DEFALT(9) 

IF(IWMEC(5,NC).NE.0) X(8)=DESVAL(8,NC) 

IF(IWMEC(5,NC).NE.O) X(9)=0ESVAL(9,NC) 

GBKL*X{8) 

GBK0=X(9) 
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C —CALCULATE G.B. DIMENSIONS — 

GBLN=0. 005*46. *(SHP**.5)*((1/S6R)**.33)*GBKL ORIQUMAL PAGE IS 

GBDI=0. 005*25. *(SHP**.5)*((1/S6R)**.33)*GBKD OF POOR QUALITY 

ILENG(NC)=GBLN 

C — ABOVE CALCS. BASED ON 1980 DDA "STAT" REPORT BASELINE ENG — 

101 CONTINUE 

WRITE( 10. 1 20)T0RK,GBWT,GBLN,GBDI 
C 
C 

C — TEST FOR CHOICE TO WEIGH PROP — 

IF(IWMEC(2,ho).EQ.O) WRITE( 10,230) 

IF(IWMEC(2,NC).EQ.O) GO TO 300 
C 

c ******** begin calcs of prop WT, ******************* 

C — CHECK FOR INPUT VALUES OF PROP WT. EQ. EXPONENTS— - 
DO 20 1=10,15 
X(I)=DEFALT(I) 

IF(IWMECI2,NC).NE.0.AND.DESVAL(2,NC).NE.0.0) X( I )=DESVAL( I,NC) 

20 CONTINUE’ 

C — DEFINE COEFFICIENTS OF PROP WT EQ — 

D=SQRT ( DATOUM ( 1 , NC ) /DATOUM( 3 , NC ) ) 

B=DATINP(5,NC) 

AF=DATINP(7,NC) 

PN=DArOUM(2,NC) 

SHP=DAT0UM(1,NC) 

PM=DESVAL(1.NC) 

WRITE(10,130)D,B,AF,PN 
IF(DESVAL(2,NC).EQ.0.0) PKW=DEFALT(2) 

I F ( DES VAL ( 2 , NC ) . NE .0 . 0 ) PKW=DESVAL ( 2 , NC ) 

C — CHECK FOR COUNTERWEIGHTS — 

IF(IWMEC(3,NC).NE.O) GO TO 200 
C 
C 

C *** BEGIN CALC OF PROP WT. (NO COUNTERWEIGHTS) ***** 

C 

C— - NOTE:PROP WT.EQ. IS BROKEN INTO THREE PARTS FOR CLARITY 

PART1=(D/10.)**X(10)*(B/4.0)**X(11)*(AF/100.)**X(12) 

PART2=((PN*D)/20000.)**X(13)*(PM+1.0)**X(14) 

PART3=(SHP/(10.*D**2))**X(15) 

PRPWT=PKW*(PART1*PART2*PART3) 

C — COMBINE PROP WT. W/ G.B. WT. — - 
WATE(NC)=PRPWT+GBWT 
WRITE(10,140)PRPWT,WATE(NC) 

C— SKIP OVER NEXT PROP WT. CALC — 

GO TO 300 
200 CONTINUE 
C 

C **** begin PROP WT. CALC FOR SINGLE ACTING PROP (W/COUNTERWEIGHTS) 

C 

CW=0.0 

C — SET EXPONENT VALUES FOR COUNTERWEIGHT WT. EQ. — 

X(12)=0.7 

X(13)=0.4 
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C — CHECK FOR INPUT WT. EQ. SCALAR{''KW VALUE) — 

C NOTE: PROP WT.EQ. IS BROKEN INTO THREE PARTS FOR CLARITY 

PART1=(D/10. )**X(10)*(B/4.0)**X(11 )*(AF/100. )**X(12) 

PART2*((PN*D)/20000.)**X(13>*(PM+1.0)**X(14) 

PART3=(SHP/(10.*l}**2))**X(lf.) 

PRPWTC=PKW*(PART1*PART2*PART3) 

C — CHECK FOR INPUT OF CON*'T. VALUE CW WT. — 
IF(DESVAL(7,NC).NE.O.O) X(/)=DESVAL(7,NC) 
IF(DESVAL(7,NC).£Q.O.O) - 

I X ( 7 ) =5* ( ( 0/ 1 0 . ) **2*B* ( AF /1 00 . )**2* ( 20000 . / ( PN*D ) ) ** . 3 ) 
CWT=X(7) 

C — SUM PROP + CW WIS. 

PRPWT=PRPWTC+CWT 

WATE(NC)=PRPWT+GBWT 

WRITE(IC,150)PRPWT,CWT,WATE(NC) 

C 

300 CONTINUE 

C 

C 

no F0RMAT('0',2X,' SHAFT-PROP GEAR RATIO®' ,F10.5,5X,- 
TMAX RPM=',F15.4,5X,'SHAFT H.P.®' ,F15.4) 

120 FORMAT( 'O', 2X, 'TORQUE®', F15.4,11X, 'GEAR BOX WT.®' .F15.4/2X, 

I'GEAR BOX LENGTH®', F15. 4, 5X,' GEAR BOX DIA.®' ,F15.4) 

130 F0RMAT('0',2X,'PR0P DIA.®' ,F15. 4, 9X, 'NUMBER BLADES®', - 

1F15.4/'ACTIVITY FACTOR®' , FI 5.4, 6X, 'PROP RPM®',F15.4) 

140 F0RMAT('0',2X,' PROP WT.®' ,F15.4,9X, 'PROP + GB WT.®',F15.4) 

150 F0RMAT('0',2X,' PROP WT.®' ,F15.4,5X, 'COUNTERWEIGHT WT.®' , - 

1F15 4 ?X 'PROP + rw WT =' F15 41 

210 F0RMAT(3X,'*** GEARBOX*WEIGHT*WILL NOT BE CALCULATED ***' ) 

220 FORMAT (3X,'*** GEARBOX DIMENSIONS WILL NOT BE CALCULATED ***') 
230 F0RMAT(3X,'*** PROPELLER WEIGHT WILL NOT BE CALCULATED ***') 

C 

IDID(NC)=1.0 

C 

RETURN 

END 


C SUBROUTINE WTEST 

C 

c 

C PURPOSE 

C 

C TO CONTROL THE CALLING OF SUBROUTINES WHICH WILL ESTIMATE THE 

C WEIGHT AND LENGTH OF INDIVIDUAL COMPONENTS 

C 

C DESCRIPTION 

C 

C THE OVERALL LENGTH OF THE ENGINE IS CALCULATED BY PROCESSING THE 
C ILENG ARRAY.ALL COMPONENTS EXCEPT DUCTS AND SHAFTS, THEN DUCTS. 

C THE REMAINING COMPONENTS EXCEPT DUCTS AND SHAFTS ARE PROCESSED. 

C THE DUCTS ARE PROCESSED AND FINIALY THE SHAFTS. 

C A BUILT-IN ASSUMPTION IN THE DUCT ROUTINE IS THAT NO DUCT IS 


25 



c -^(— (Tj r~) r~>c~> c~5r~>nc~>r~jc^<n nnr^r^c^<r> r^<~>c~vc'^r^r^r^ r^<~><~5r~>r^r~ior~>r~>c~> 


CONNECTED TO ANOTHER DUCT I.E. THE DUCT SIZE IS DETERMINED BY THE 
ADJOINING COMPONENTS. 

THEN THE MAXIMUM RADIUS IS FOUND. THEN DEPENDING ON THE PRINT 
FLAG -lO'JTCD- THE REQUIRED PRINTING IS DONE 
IT THE PLOT CODE FLAG -IPLT- IS TRUE ROUTINE EPLT IS CALLED 
USAGE 


CALL WTEST 
CALLING ROUTINES 


FLOCAL- 
ZTOPZ - 

REQUIRED SUBROUTINES ORIGiNAL PAGE ?S 

OF POOR QUALITY 

COMP -COMPRESSOR WEIGHT/LENGTH 
TURB -TURBINE 
SHAFT -SHAFT 
DUCTW -DUCT 

COMBWT -PRIMARY BURNER WEIGHT/LENGTH 
WTNOZ -NOZZLE WEIGHT/LENGTH 
WMIXR -MIXER 
WSPLT -SPLITTER 
EPLT -PRINTER/PLOTTER 

MODIFICATION HISTORY 


DATE ID ANALYST 


DESCRIPTION 


MO/DA/YR IDENT NAME DESCRIPTION OF CHANGES 

9/21/Sl ROBERT COPBAN INCORPORATE INTERACTIVE MODE AND GR 

Hies OPTION 

8/21/82 PETE SENTY INCORPORATE PROPELLER WEIGHT MODEL 

AUTHOR/LANGUAGE/DATE 


NORMAN PREWITT-BOEING COMPUTER SERV. /FORTRAN IV / OCT 10,1976 
GLOSSARY 


NAME ORIGIN USAGE DESCRIPTION 


C 


IWMEC 

/WMECH/ 

I 

WATE 

/WMECH/ 

0 

ALENG 

/WMECH/ 

0 

TLENG 

/WMECH/ 

0 

RI 

/WMECH/ 

0 

RO 

/WMECH/ 

0 

DESVAL 

/WMECH/ 

I 

DSHAF 

/WMECH/ 

0 

RPMT 

/WMECH/ 

I 

IWT 

/WMECH/ 

I 


CONTROL INFORMATION 
WEIGHT OF EACH COMPONENT 
ACTUAL LENGTH OF EACH COMPONENT 
ACCUMULATED LENGTH TO END OF COMPONENT 
RADIUS INNER INLET, OUTLET EACH STATION 
RADIUS OUTER INLET, OUTLET EACH STATION 
MECHANICAL DESIGN DATA OVERRIDES DEFAUL 
SHAFT DIAMETER INNER TO OUTER 
ACTUAL COMPONENT RPM 
WEIGHT ESTIMATION FLAG TRU= DO IT 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 


c 

c 

c 


IPLT /WMECH/ I 
PLOT L 
lERR /WMECH/ 0 

1510 /WMECH/ 1 

1511 /WMECH/ I 
lOUTCD /WMECH/ I 
ILENG /WMECH/ I 
IDID 0 


f^OOR 

PLOTTER FLAG TRUE= DO IT 
GRAPHICS PLOTTING FLAG TRUE* DO IT 
ERROR FLAG 

OUTPUT UNITS 0=-EN6LISH, 0 SI 
INPUT UNITS 0=EN6LISH, 0 SI 



PRINT FLAG 0=SUMMARY, l=GENERAL.2»DIACNOSIIC 
COMPONENTS CONTRIBUTING TO OVERALL LENGTH 


FLAG = 0 COMPONENT NOT YET WEIGHED =1 YES 


SKI PIT L BACKGROUND FLAG YES^SKIP PROMPTS 

**1t**1t1i1t*1r**1t***K*1t**1t****1i*1t*if1rk***1rk**1l1Ht******1HHt*ik**1t1t1t*1Ht*t** 


SUBROUTINE WTEST 

REAL*8 DAT I NP , DATOUT . WTF , TOPRES , TOTEMP , FAR , CORFLO , VMACH . ST ATP , ERRO- 
1 R, TOL , TOLT , TOLTT . DEPV , DTOL , PERPF .RPMNT, TMTEMP , TMPRES , DATCUM, UATMAC - 

2. DATALT,DESLIM,TNPRES,TNTEMP,CNRFL0.C0RFLM,WNTF,DAT0UN,FARN,DANINP- 

3 , DEBUG , DEPQ , SELAST , DD , TOLOPT , RCH . STOC , TFUEL 

***************** 

* COMMON BLOCKS * 

***************** 

COMMON /IJBL/ DAT I NP ( 1 5 . 60 ) . DATOUT ( 9 , 60 ) , WTF ( 40 ) , TOPRES ( 40 ) , TOTEMP ( - 
140),FAR(40),C0RFL0(40).VMACH(40),STATP(40),ERR0R(40),T0L,T0LT,T0LT- 
F T , DEPV ( FO ), D TOL ( 20 ), PERPF ( 20 ), RCH . STOC , TFUEL 
COMMON /SNGL/JM1,JM2,JPI,JP2,JCX,L0CTBL(9,60),JC0MP(70),IWAY.NIT,- 
IITAB(/0),JC0NF(6C,4;,JTYPE(60),JFL0W(70),IDEDAP(15).KKINDS(14,25),- 
2NC0MP,N0STAT,NITER,NFINIS.NPASS,JCC,NrBL,NCTS,JCIND(20),JCDEP(20).- 
3JCVIND(20),JCVDEP(20).KDTYP(20),ID0NE(60) 

COMMON /DEPAUL/ DEFAUL( 15,20) , 1SCALE(3),SCALE(6) 

COMMON /WMECH/ IWMEC(7,60),WATE(60),ALENG(60),TLENG(40).RI(2,40),R- 
i0{2, 40), UESVAL( 17,60), DSHAF(5),RPMT(60),IWT,IPLT,IERR, ISII, ISIO, 10- 
2UTCU.NSrAG(60) 

COMMON /CONVER/ C0NVFR(15) 

COMMON /NEPOPT/ DEBUG , DEPQ . SELAST , DD , TOLOPT , NDSET, NPARTS , lOPTP , NPA- 
ISSO,NVOPT,NJOPT,NOPT 

COMMON /TERMON/ TNPRES ( 40 ) , TNTEMP ( 40 ) , CNRFLO ( 40 ) , WNTF ( 40 ) , RPMNT ( 40- 
1 ) . DATOUN (9,60), FARN ( 40 ) , DAN I NP ( 1 5 ,60 ) , TMTEMP ( 40 ) , TMPRES(40 ) , CORFLM- 
2{40),DAT0UM(9,60),DATMAC(4,60),DATALT(4,60),DESLIM{15) 

COMMON /CENTER/ CGARM(60) 

COMMON /ARM/ ACCARM 
COMMON /DISKK/ DISKWI , ENGINE 
COMMON /SKIP/ SKIPIT 
COMMON /WTLDAT/ WATENG,WATACC 

* DATA STORAGE DEFINITION * 

LOGICAL PINP,IPLT, ISIO, ISII, PLOT, SKIPIT 
INTEGER IDID(60),ILENG(40) 

DIMENSION NUMNUM(17), IRNAME(17), C0RFLC(40) 

NAMELIST /W/ IWMEC.DESVAL,ACCS,IWT,IPLT,ISII,ISI0,I0UTCD,ILEN6,DES- 
1LIM,ISCALE, SCALE, ACCARM, DISKWI, ENGINE, PLOT, IDID 

* DAIA STATEMENTS * 

'kik'kitif'kii'kA ' \ -k it 'k iie it it if 
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DAI A UUII-.ISIIAI .hNliU,SUI.lVALV/‘UUHU:i,4HSHAF.4!itNUL..4HSUl .4HVALV/ 
DA I A 1 U'l‘ . I HIH: , 11 AN . U 0 . 1 n , 1 HD r , 1 1 r I / JML , .DIMIH: , JHl AN . .*HFO . ,‘Hl 1 , ■ 

MllHt‘1 ,.UIl Dl/ 

DAIA IMNI‘.YlS/.IKDt..4HY /.DLOI/.FALSt./ 

ItSI Wits I KLAU 

IF (IWI.IQ.D) HO 10 tn‘0 
IF IIWAY.HE.O) HO 10 H) 

DO 40 1-1,40 
WNIFlD-WMlD 
FAKNtD-l AKll) 

INl’KlSlD-IODKtSlD 

inii.mi*(1)-ioii:mi’(1) 

lNKFl.O(l)-a)RF 1.00) 

DA100N(1.1)-DAUUI1(I,1) 

DAlOUNlD.n-DAlOUIU'.l) 

DAIOUNIO.D-DAIOUUH.I) 

DA 1 00N( 4,1) -DA 100 1(4,1) 

U (l.U./O) HO ID 10 
K-1*. 0 

DAI00N(1,K)-DA10UU1,K 
DA100N(. ,K)-DA100T(,’,K 
l'A10DN({5,k)-DA1001(0,K) 

DA100N'4,R)-DAIOD1(4,R) 

10 IMDKt.S(l)H) 

INllMI'i D-O 
lOKllM( l)-0 
DO A) KK-1,4 
DAIMAi'(KK,l)-0 
A) DAlAl l(KK,l)-0 
DO .50 KK-1 ,4 
DA10DM(KK,l)-0 
It (1.1 I ..’O) HO 10 .to 
K- 1*.’0 

DA10DM(KK,K)-0 
.to I'ONIINOI 
40 l ONl INOl. 

DO !h) 1-1,400 
:>0 DAN1NI’(1,1)-DA11NI'(1,1) 
t»0 DO ttO ;-l,40 

4 A A A / • . * 

It (OlYl'l (!).N1 .10) HO 10 
It (DAI lNI’(.‘,l).i:g.O.O) HO 10 

It I DA 1 ODM ( 1 , 1 ) . I 1 . D Alt S ( DA 1 00 1 ( 1 , 1 ) ) ) DA lODM ( .t , I ) — D A 1 (HI I ( 1 , 1 ) * 4 . {1 • 
1 t)4oH * DA 1 00 1 ( .' , I ) * * .V ( .ttHtO . * DA 1 00 1 ( / , I ) * V ) 

I i (DAI DOM ( I , I ) . I I . DAlt S ( DA 1001(1.1))) DAI OOM ( . 1 ) -DAI 00 1 ( .',1 ) 

^ 1 1 ( DA 1 OOM ( 1 . 1 ) . I 1 . DAlt S ( da 1 00 1(1,1))) DA 1 OOM ( 1 . 1 ) -DAWS ( DA 1 00 1(1,1). 

ot> I'ONIINOL 

4 444 4 

I'OKt U:( 1)-1'0KF10( l)/l .S4‘'/.'M> 

It (.llYt'l (D.llt.l) ll-l 


.\'t 


k 



II- (roi’KtS(i) 
IF (rOPKtS(l) 
IF (TOTtMPII 
IF (TOTEMP(.I) 
IF (UATUinic, 
IF {l)ATOUT(.\ 
IF (i:ORFl.i:(l) 
IF (CORFLC(l) 
IF irOPRES(l) 
IF (rOTEMP(l) 
IF (CORFLC(l) 
IF (DATOUT(P, 


.GT.lMPRtS(I)) 

.GT.TMPRES(l)) 

.GT.lMTEMPIl 

.GT.TMrEMP(l)) 

n.6T.l)AT0UM(:.’ 

l).GT.DATOUM(^ 

.GT.CORFLM(l)) 

.GT.CORFLM(l)) 

.GT.TMPRESil 

.GT.TMTEMP(l)) 

.GT.CORFLM(l)) 

I).GT.UATOUM(:’ 


ORlGmi 

OP POOR 


UArMAC(l,l)»DATlNP(b,II) 
DATALr(l,l)»DATlNP(9.1lj 
l)ATMAC(?.I)-0ATINP(5,n) 
DArALl(i?.I)=l)AriNP(9.1I) 
,1)) DATMAC(3.I)-0ATlNP(b,II 
.1)) DATALT(3.1)-DAT1NP(9,I1 
l)ATMAC(4,I)»l)ATlNP(b,II) 
t)ATALT(4,I)-DATINP(9.II) 
TMPRESII »TOPRES(l) 
TMTEMP(l)»TOTEMP(n 
(;ORFLM(I)=CORFLf(I) 

,1)) DATOUI(2,1)»DA10UT(^,1) 


IF (l.LE.20) GO ro 70 

IF (0AT0UT(2,K).GT.DAT0UM(2,K)) 0Ar0UM(2.K)-l)AT0UT(2,K) 

70 IF (JTVPE(1).NE.5) GO TO 80 
UATTRQ=OATOUT( I , I )/DAT0UT(2, I ) 

IF (UATTRQ.GT.DATOUMIl.l)) DATOUMI 3. 1)=0AT0UT(2, I ) 

IF (i)ATTKQ.GT.nATOUM(lJ)) l)ATOUM( 1 , 1 )=OATTRQ 
IF (l.LE.20.0R.JTYPE(K).NE.5) GO TO 80 
l)ATTRQ=U)ATOUT( 1 .K)/8AT0UT(2.K) 

IF (1)ATTRI}.GT.DAT01IM(1.K)) l)ATOUM(3,K)“OATOUT(2,K ) 

IF (l)ATTRQ.GT.0A10UM(l,K)) DATOUM( I ,K)»8ArTRQ 

THE FOLLOWING MAY BE USEFUL FOR DEBUGGING AND HAS BEEN LEFT 1 

WRITE(IO,9y2)jrYPE(l) 

992 FORMAT (//' COMP TYPE *13) 

WRI1E( 10,997) 

99/ FORMATC DES lOIPRES N 0 rOFEMP') 

WR 1 1 E 1 1 0 , 999 ) TORRES ( I ) , TNPRE S ( 1 ) . 1 0 TEMP (I ) . TN TEMP ( I ) 

WRITE (10,998) 

99t> F'ORMATC DES CORFLO N D M WTF') 

WRITE! 10, 999 ) CORF L0( I ),CNRFlO( 1 ),C0RFLM( I ),WNTF( 1 ) 

WR1IE( 10,993) 

993 FORMAU* DES EFF N D PR') 

WRl rE( 10,999)DAI0UT(8, 1 ) ,DA10UN(8. 1 ) ,DAT0UT(9, 1 ) ,DAT0UN(9, 1 ) 
WRIIEl 10,99b) 

99'j FORMAT!' MAX TO I PRES N M TOTEMP') 

WRl IE( 10,999) 10PRES( 1 ),TMPRES( 1 ), TOTEMP ( 1 ),TMTEMP(I ) 

WRl IE 1 10,994) 

994 FORMATC MAX RPM N M D') 

WR 1 TL ( 1 0 , 999 ) DA I Oil T ( .' , 1 ) , DATOUM( 2 , 1 ) , DAT OUN ( 2 , 1 ) 

WRITE! 10,991) 

991 FORMAU ' MTRQ N M TRPM MIRPM ') 

WR 1 TL ( 1 0 , 999 ) DA T 1 RQ , DA TOUM (1,1), DATOU T ( 2 , 1 ) , DATOUM( 3,1) 

999 FORMAT (4F 12.3) 

kM 1 1 E ( 10 , 999 ) DATMAC (1,1), DATMAC ( 2 , 1 ) . DA TMAC (3.1) 

WRl TE( 10,999)DA1ALT( 1 , I ),DATALT(2, 1 ) ,DATALT(3, I ) 

80 C0N11NUE 

IF (IWT.GE.2) GO TO 90 

ZERO OUT OUTPUT ARRAYS 

GO TO li20 
90 JSCALE=0 
1SAVE=1W1 

IF (1WT.NE.4) GO TO 130 


WR 1 TL ( 1 0 , 999 ) DA I Oil T ( 2 . 1 ) . DATOUM( 2 . 1 ) . DAT OUN ( 2 . 1 ) 



C-) <-> C~> O 


100 JSCALE=JSCALE+1 

IF (JSCALE.GT.ISCALE(2)) 60 TO 620 

IF (JSCALE.GT.1.AND.ISCALE(1).EQ.2) I0UTCD=1 

SCALEF=SCALE(JSCALE) 

IF (JSCALE.GT.l) SCALEF=SCALE( JSCALE)/SCALE( JSCALE-1 ) 
no IF (REVISE. EQ. YES) SCALEF=1 .0 
DO 120 1=1,40 
WNTF(I)=WNTF(I)*SCALEF 
CNRFLO( I )=CNRFLO( I )*SCALEF 
CORFLM( I )=CORFLM( I )*SCALEF 


IDID(I)=0 

DATOUN { 1 , 1 ) =DATOUN (1,1 )*SCALEF 
CONTINUE 

IF (JSCALE.GT.l. AND. IWT.EQ. 4) GO TO 190 
IF (REVISE. EQ. YES) 60 TO 180 
130 00 140 1=1,5 
140 DSHAF(I)=0. 

DO 150 1=1,60 

WATE(I)=0 

NSTAG(I)=0 

IDID(I)=0 

RPMT(I)=0. 

150 ALENG(I)=0 


120 
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DO 160 1=1,40 
TLENG(I)=0 
RI(1,I)=0 
ILENG(I)=0 
RI(2,I)=0 
R0(1,I)=0 
160 R0(2,I)=0 


NAMELIST READ OF WTEST DATA 


CALL NAMEPR (9,10,8,PINP) 

READ (8,W) 

DO 170 1=1,60 
DO 170 K=l,17 

IF (IWMEC(1,I).EQ.NUMNUM(K)) IWMEC( 1 , I )=IRNAME(K) 

170 CONTINUE 
GO TO 190 
180 WRITE (20,760) 

CALL NAMEPR (20,10,8,PINP) 

READ (8,W) 

PROCESS LENGTH CONTRIBUTING VECTOR EXCEPT DUCTS AND SHAFTS 

190 DO 310 1=1,40 
NC=ILENG(I) 

C WRITE(10,7777)NC,ILEN6(I),JTYPE(NC) 

IF (NC.EQ.O) GO TO 320 
JT=JTYPE(NC) 

60 TO (310,230,270,200.210,280,260,220,250,295,310.310,310,310). J- 
IT 

C COMPRESSOR 


30 



200 CALL COMP (NC) 

GO TO 300 

C TURBINE ORIGJNAL PAGE \S 

210 CALL TURB (NC) OF POOR QUALITY 

GO TO 300 

C MIXER 

220 CALL WMIXR (NC) 

GO TO 300 

C PRIMARY BURNER 

^30 IF (IWMEC(l.NC).EQ.IDUC) GO TO 240 
IF (IWMEC(l.NC).EQ.IVALV) GO TO 290 
CALL COMBWT (NC) 

GO TO 300 

C — DUCTS 
240 CALL DUCTW (NC) 

GO TO 300 

C NOZZLES 

2*^0 CALL WTNOZ (NC) 

GO TO 300 

C SPLITTER 

260 CALL WSPLT (NC) 

GO TO 300 

C TRANSFER DIMENSIONS FOR WATER INJECTION 
270 CALL DUMMY (NC) 

GO TO 300 

C HEAT EXCHANGER WEIGHT 
280 CALL HMEC (NC) 

GO TO 300 

C VALVES 

290 CALL VALVWT (NC) 

********************** 07/15/82 ***** 

GO TO 300 

295 IF (DATINP(2,NC).EQ.O.O) GO TO 310 
CALL PROPWT (NC) 

IUP=JC0NF(1,1) 

1DN=JC0NF(1,3) 

GO TO 305 

0 ********************** 07/15/82 ***** 

C ACCUME LENGTH 

300 IUP=JCONF(NC.l) 

IDN=JC0NF(NC,3) 

305 TLENG( I DN ) =TLENG( I UP )+ALENG( NC ) 

IF (JT.EQ.6) TLENG(IDN)=TLENG(IUP) 

ID2=JC0NF(NC,4) 

IF (ID2.GT.0) TLENG(ID2)=TLENG(IUP)+ALEN6(NC) 

IF (JT.EQ.6) TLENG(ID2)=TLENG(IUP) 

IU2=JC0NF(NC,2) 

IF (IU2.GT.0) TLENG(IU2)=TLENG(IUP) 

IDID(NC)=1 
310 CONTINUE 

C LAST COMPONENT WAS A NOZZLE SET ENGINE MAXIMUM LENGTH 

320 ENGLEN=TLENG(IDN) 

C 
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C PROCESS REMAINING COMPONENTS 

C 

00 410 1=1,60 

C WRITE(10,7777)NC,1,IDI0(I) 

IF (lOlD(I).EQ.l) 60 TO 410 
C 

C - PROCESS COMPRESSORS, TURBINES, MIXERS, BURNERS, SPLITTERS 
NC=JTYPE(I) 

IF (NC.LE.O) 60 TO 410 

GO TO (410,360,400,330,340,390,380.350,410.410,410,410,410.410), N- 

1 c 

C — COMPRESSOR 
330 CALL COMP (I) 

GO TO 410 

C — TURBINES 
340 CALL TURB (I) 

GO TO 410 
C — MIXER 
350 CALL WMIXR (I) 

GO TO 410 

C — BURNERS 

360 IF (IWMEC{I,I).EQ.I0UC) GO TO 110 
IF (IWMEC(1,I).EQ.1VALV) GO 10 370 
CALL COMBWT (I) 

GO TO 410 

I' VAI UFS 

370 CALL VALVWT (I) 

GO TO 410 

C --- SPLITTER 
380 CALL WSPLT (I) 

GO TO 410 

C HEAT EXCHANGERS 
390 CALL HMEC (I) 

GO TO 410 

C TRANSFER DIMENSIONS FOR WATER INJECTION 
400 CALL DUMMY (I) 

410 CONTINUE 

C WR1TE(10,7771)JTYPE 
C 

C PROCESS DUCTS 

DO 4^0 1=1,60 

IF (IDID(I).EQ.l) GO TO 420 
NC=JTYPE(I) 

IF (NC.NE.2) GO TO 420 
IF (IWMEC(1,I).NE.IDUC) GO TO 420 
CALL DUCTW (I) 

420 CONTINUE 

C WRITE(10,7771)JTYPE 

C PROCESS NOZZLES 

DO 430 1=1,60 

IF (^IDID(I).EQ.l) GO TO 430 
NC=JTYPE(I) 

IF (NC.NE.9) GO TO 430 
CALL WTNOZ (I) 

430 CONTINUE 
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C WRITE(10,7?7UOTYPE ORIGINAL PAGE IS 

C - — ACCUME itNGTH OF POOR QUALITY 

00 450 t*\,40 

C WRITE()0«7777)NC,I,tDID(I) 

IF (JTVPE(l).LE.O) GO TO 440 
IF (lOIO(I).EQJ) GO TO 450 
NC-OTYPE(I) 

GO TO (450.440*440, 440, 440, 440, 440, 440, 440, 450,450, 450,450, 450) » N- 
IC 

440 IUP«JCONF(I,n 
IU2-OCONF(I,2) 

IDN-JC0NF(I,3) 

ID2-JCONF(t,4) 

TLENG( IDN)-TLENG( IUP)’^ALENG( I) 

IF (NC,EQ.6) TLENG(I0N)«TLENG(IUP) 

IF (IU2.GT.6) TLENG(1U2)«TLENG(IUP) 

IF (1D2,GT.0) TIENG(I02}«TLENG(I0N) 

IDID(I)-1 
450 CONTINUE 

C PROCESS SHAFTS 

DO 470 >1,5 
DO 460 1*1,25 
NC*KKINDS(U,I) 

IF (NC.LE.O) GO TO 470 
IF (INNEC(\,NC).NE.LSHAF) GO TO 460 
IF (IWNEC(2,NC).EQ,J) CALL SHAFT (NC) 

460 CONTINUE 
470 CONTINUE 
C 

C - — FIND ENGINE NAXIMUM RADIUS 

XR-0 

DO 490 I-I,NOSTAT 

IF (XR.GE.RO(M)) GO TO 480 

XR-RO(M) 

480 IF (XR.GE.K0(2,t)) GO TO 490 
XR*R0(2.I) 

490 CONTINUE 

C 

C • — GET ENGINE TOTAL WEIGHT AND ALEN6 CONVERSION 

WATENG»0 i 

IF (ACCS.EQ.O) ACCS*.! 

HAT-0. 

DO 500 1-1,60 

IF (vlTYPE(I).EQ.9) GO TO 500 
WAT*WATE{I)+WAT 
SOO CONTINUE 

WATACC-ACCS^WAT 

IF (lOUTCD.taU) WRITE (10,660) lOUTCD 
IF (ISIO) HArACC*WATACC*C0NVER(3) 

IF (lOUTCD.Gl.l) WRITE (10,690) WATACC 
DO 510 1*1,60 
WFACTH"! ft 

IF (0ESVAL(15.I).NE,0.) WFACTR*0ESVAL(15.I) 

HATE(1)*WATE(I)*HFACTR 


^ 1 
li. i, 
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c 

c 


c 

c 


c 


IF ((JTYPE(I).EQ.4.0R.JTYPE(I).EQ.5).AND.WATE(I).EQ.0.) WATE(I)*+1- 
1.E6 

IF (.NOT.ISIO) 60 TO 510 
WATE(n=WATE(I)*C0NVER(3) 

ALER6(i)=A!.ENG(I'*C0NVER(l) 

5 1 0 WATENG=WATEN6+WATE ( I ) 


CONVERT RADIAL DIMENSIONS AND TLENG 

IF (.NOT.ISIO) GO TO 530 

DO 520 I=1,N0STAT 

RI(1,I)=RI(1,I)*C0NVER(1) 

RI(2,I)=RI(2.I)*C0NVER(1) 

R0(1,I)=R0(1,I)*C0NVER(1) 

R0(2,I)=R0(2,I)*C0NVER(1) 

TLENG( I )=TLENG( I )*CONVER( 1 ) 

520 CONTINUE 


ORIGINAL PAGE ff? 
OF POOR OUALiTv 


WRITE COMPONENT WEIGHT INFO 

530 UNITSI=ENGU 

IF (ISII) UNITSI=SIU 
UNITSO=ENGU 
IF (ISIO) UNITSO=SIU 
WRITE (10,670) UNITSI,UNITSO 
IF (IWT.EQ.4) WRITE (10,630) 

IF (IWT.EQ.4) WRITE (10,640) 

IF (IWT.EQ.4) WRITE (10,650) SCALE( JSCALE) 

IF (IWT.EQ.4) WRITE (10,640) 

IF (IOUTCD.lt. 1) 60 TO 580 
WRITE (10,680) 

HLENG=0. 

CGLEN6=0. 

C6WATE=0. 

CGT0TM=0. 

CGC0MP=0. 

DO 570 1=1,60 
NC=JTYPE(I) 

IF (NC.LE.O) GO TO 570 
WRITE(10,7777)JT 

GO TO (540,540,540,540,540,540,540,540,540,570,540,570,570,570), N- 
IC 


0 **************** 07/15/82 ********* 

540 IF (NC.EQ.10.ANO.OATINP(2,I).EQ.O.O) GO TO 570 
IUP1=JC0NF(I,1) 

IUP2=JCONF(I,2) 

IDN1=JC0NF(1,3) 

I0N2=JC0NF(I,4) 

IF (NC.EQ.IO) IUP1=JC0NF(I,1) 

IF (NC.EQ.IO) IUP2=OCONF(I,2) 

IF (NC.EQ.IO) IDN1=JC0NF(I,3) 

IF (NC.EQ.IO) IDN2=JC0NF(I,4) 

********** 07/15/82 ****** 

WRITE (10,700) I,WATE(I),ALEN6(I),TLEN6(I0N1),RI(1,IUP1),R0(1,IUP1- 

1) ,RI(1,IUP2),R0(1,IUP2),RI(2,IDN1),R0(2,I0N1),RI(2,IDN2),R0(2,IDN2- 

2) ,NSTAG(I) 
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OF POOR 0-''5-'TY 


IF (ALENGd).EQ.O.) GO TO 550 
HLENG=ALENG(I)/2. 

IF (IWMECil,I).EQ.ILPC.OR.IWMEC(l,I).EQ.IHPC) CGX»ALENG(I)-CGARM(I- 

^IF (IWMEC(1,I).EQ.IFAN.0R.IWMEC(1,I).EQ.IF0) CGX=ALENG(I)-CGARM(I) 
IF ( IWMFC( 1 , I ) .EQ. IFI .OR. IWMEC( 1 , I ) .EQ. IHPT) CGX=ALENG( I )-CGARM{ I ) 
IF (lwI*iEC(l,I).EQ.ILPT) CGX=ALENG(I)-CGARM(I) 

IF (IWMEC(1,I).EQ.ILPC.0R.IWMEC(1,I).EQ.IHPC) HLENG=CGX 
IF (IWMEC 1,I).EQ.IFAN.0R.IWMEC(1,I).EQ.IF0) HLENG=CGX 
IF (IWMEC(1,I).EQ.IFI.0R.IWMEC(1.I).EQ.IHPT) HLENG*CGX 
IF (IWMEC(1,I).EQ.ILPT) HLENG=CGX 
CGLENG=TLENG(IDN1)-HLENG 
GO TO 560 

550 IF (IWMEC(l.I).EQ.ISHAF) CGLENG=CGARM(I) 

IF (IWMEC(l,lj.NE.ISHAF) GO TO 570 

560 CGWATE=WATE(I) 

CGCOMP=CGWATE*CGLENG 

CGTOTM=CGTOTM+CGCOMP 


570 CONTINUE 

CENGR A= ( CGTOTM+ ( WATACC*ACCARM ) ) /WATENG 

— make summary print 

580 IF (.NOT.ISIO) GO TO 590 
XR=XR*C0NVER(1) 

ENGLEN=ENGLEN*C0NVER(1) 

CENGRA=CENGRA*C0NVER(1 ) 

590 WRITE (10,710) WATENG, WATACC.ENGLEN.XR 
IF ( SCALE (JSCALE).EQ.l.) SEXP01=WATENG 
SEXP02=1 . 

IF ( SCALE ( JSCALE ) . NE . 1 . . AND . I WT . EQ . 4 ) SE XPOE=ALOG ( WATENG/SEXPOl ) /A- 
1 LOG ( SCALE ( JSCALE ) /SEXP02 ) 

IF ( SCALE ( JSCALE ).EQ.l.) SEXP0E=1. 

WRITE (10,720) CENGRA 

IF (IWT.EQ.4) WRITE (10,730) SEXPOE 

IF (JSCALE. GT.l. AND. ISCALE(1).EQ. 2) GO TO 610 

IF (ENGINE. EQ. 2.) GO TO 610 

IF (IPLT) CALL ENGPLT (ENGLEN.XR) 

IF (PLOT) GO TO 600 
GO TO 610 

600 IF (SKI PIT) CALL EGPLOT (ENGLEN, WATENG) 

IF (SKIPIT) GO TO 610 
WRITE (20,740) 

READ (20,770) ANSWER 

IF ( ANSWER. NE. YES) GO TO 610 

CALL EGPLOT (ENGLEN, WATENG) 

WRITE (20,750) 

READ (20,770) REVISE 
IF (REVISE. EQ. YES) IWT*2 
610 ENGINE=1. 

IF (REVISE. EQ. YES) GO TO 110 
IWT=ISAVE 

IF (IWT.EQ.4) GO TO 100 
620 IWT=0 
I0UTCD=2 
RETURN 
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ORIGINAL PAGE IS 

^ OF POOR QUALITY 

630 FORMAT (IH /24H ENGINE SCALING DATA ) 

640 FORMAT (20H *******************) 

650 FORMAT (14H SCALE FACTOR ,F5.2) 

660 FORMAT (IH /14H *************/ 14H * */,14H * ACCS WT *- 

1/14H * */13H ************, 1 ;^ 

670 FORMAT (1H1.26H WEIGHT INPUT DATA IN ,A4,6H UNITS/27H WEIG- 
IHT OUTPUT DATA IN ,A4,6H UNITS//) 

680 FORMAT (69H COMP WT COMP ACCU UPSTREAM RADIUS DOWNS- 
1TREAM RADIUS /77H NO EST LEN LEN RI RO RI RO - 
2 RI RO RI RO NSTAGE/) 

690 FORMAT (/,11H ACCS WT=,F8.3) 

700 FORMAT (I7,F6.0,F7.0,F6.0,4F5.0,F6.0,3F5.0,I8) 

710 FORMAT (/,27H TOTAL BARE ENGINE WEIGHT=,F6.0,2X,12HACCESS0RIES=,F- 
17.2,2X,23HESTIMATED TOTAL LENGTH=,F6.0,2X,25HESTIMATED MAXIMUM RAD- 
2IUS=,F5.0) 

720 FORMAT (30H ESTIMATED CENTER OF GRAVITY=,F6.0) 

730 FORMAT (39H ESTIMATED AIRFLOW SCALING EXPONENT IS.F6.3) 

740 FORMAT (46H DO YOU WISH A GRAPHICS PICTURE? YES=Y ; NO=N) 

750 FORMAT (54H DO YOU WISH TO MAKE CHANGES TO THE INPUT? YES=Y;NO=N) 
760 FORMAT (42H ENTER DESIRED CHANGES FROM TERMINAL, I .E./39H &W DE- 
1 SVAL (1,6) =**** ( DATA ) **** &END ) 

770 FORMAT (A4) 

END 

'k'kic'k'k 
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TABLE I. - USER INPUT FOR PROPELLER PERFORMANCE 


SPEC 

Option 1, 
fixed efficiency 

Option 2, 
generalized map 

Option 3, 
user-input map 

1 

(-) Shaft horsepower 

(-) Shaft horsepower 

(-) Shaft horsepower 

2 

Efficiency 

Any negative number 

Map number 

3 

Static T/SHP 

— 

— 

4 

Power loading^ 

Power loading 

Power loading 

5 

Number of blades^ 

Number of blades 

Number of blades® 

6 

— 

Lift coefficient 

— 

7 

Activity factor^ 

Activity factor 

Activity factor^ 

8 

— 

Optional scaling 
factor 

Optional scaling 
factor 

9 

— 

Optional desired 
design efficiency 

Optional desired 
design efficiency 

10 

Tip speed^ 

Tip speed 

Tip speed 


aRequired only for weight calculation, not for performance. 


ORIGINAL PAGE IS 
OF POOR QUALITY 


38 




ORIGfNAL PAGE IS 
OF POOR QUALITY 



Figure L - Baseline gearbox weight as a function of 
torque for gear ratio - ^ il8. 



